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Abstract 



C*2 

£>y In the mechanics of inviscid conservative fluids, it is classical to generate the equa- 

tions of dynamics by formulating with adequate variables, that the pressure integral 
calculated in the time-space domain corresponding to the motion of the continuous 
medium is stationary. The present study extends this principle to the dynamics of 
large deformations for isentropic motions in thermo-elastic bodies: we use a new 
way of writing the equations of motion in terms of potentials and we substitute the 
trace of the stress tensor for the pressure term. 
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1 Introduction 



HAMILTON'S PRINCIPLE holds for all conservative mechanical systems with 
holonomic side conditions. This is the case with perfect fluids in adiabatic 
motion. If (F) represents the set of the virtual motions which assign to the 
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system two given positions at times t\ and £2, the principle can be stated as 
follows {principle of least action): 

Among the (F) -motions, the motion of the system is that which stationarizes 
the integral of the difference between kinetic energy and potential energy over 
the time-space domain occupied by the continuous medium. 

Several authors have observed that in conventional fluid mechanics, many 
variational principles can be reduced to Hamilton's principle, although the 
procedures are diverse, according to the choice of the unknown functions. In 
the Eulerian description of fluids, several authors [1-5] have shown that the 
equations of motion may be derived from a variational principle through the 
introduction of Lagrange multipliers corresponding to the side conditions that 
the variations of the mass density, the entropy and the Lagrangian coordinates 
are coupled with some conservation conditions. Minor details apart, and with 
V, p, s, X, 0, ip, e as independent variables, these methods consist in station- 
arizing 
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where X 



X 1 
X 2 
X 3 



are the Lagrangian coordinates, x 



x 



X 



X 



the Eulerian coor- 



dinates, p is the mass density, V the velocity, e the internal energy density, 
s the entropy density, Q the extraneous force potential; and -0 are scalar 
Lagrange multipliers, E is a vector Lagrange multiplier and T stands for the 
transposition. The variations with respect to 0,-0,3 give the constraints: 



dp 

m 

dps 
~dt 



+ divpV = 0, 



density equation 



d( P X 7 
dt 



+ divpsV = 0, entropy equation 



div(pVX r ) = 0, Lin's constraint [5,6]. 



(2) 



(3) 
(4) 



The last equation expresses that, in the eulerian description, the initial co- 
ordinates X do not change along a particle path [2]. The connection of this 
principle with Hamilton's principle is somewhat intricate for any direct use. 
Casal [7], Seliger & Whitham [8] have introduced considerable simplification. 
By integrating by parts, and supposing that the integral corresponding to the 



boundary of the time-space domain has a zero- variation, these authors reduce 
the integrand in (1) to: 

-fjV 2 — pe — pVt — pcj) — psi/j — pX r H 

where " is the material derivative. By using Clebsch representation [9] for 
the motion and thermodynamic relations, they reduce the Lagrangian density 
to be just the pressure p. Conversely they have to find a suitable variational 
principle which generates the dynamical equations of a fluid and auxiliary 
conditions by expressing that 

/ p dxdt (5) 

Jw 

is stationary. 

On this purpose, Seliger and Whitham express the pressure as a function of 
the enthalpy and the entropy. They give a Clesch representation of the velocity 
field by means of Lagrange multipliers </>, ijj, 3. This does not restrict the form 
of the velocity vector and gives a redundant number of unknowns. Under the 
constraint 

p = -pV 2 — pe — pVt — p (f) — ps ijj — /9X T S, 

the variational principle which states that (5) is stationary with respect to the 
variables <f>, s, ip, X and H , yields the equations of the fluid motion in terms 
of potentials, and the density equation. 

Obviously, this principle is far from Hamilton's principle. Seliger and Whitham 
noticed that Clesch's representation of the velocity is a crucial step for the final 
form (5) and said this seems to be an especially simple form limited to fluids. 

The purpose of the present study is to generalize this principle to isentropic 
motions in thermo-elastic bodies with thermodynamic potentials [10]. The 
medium may be inhomogeneous. 

We analyze the article by Seliger & Whitham and the limitation caused by 
Clesch representation. By using an appropriate representation for the internal 
energy, we overcome the difficulties associated with the fact that the stress 
tensor of an elastic medium is no more spherical. 

In the second section, a decomposition of the stress tensor into spherical and 
deviatoric parts is associated to partial derivations of the internal energy. The 
decomposition we obtained in 1986 is now largely used in the literature. 

In the third section, an Eulerian approach of the variational principle similar to 
the method used by Seliger and Whitham leads to a new Clesch representation 
and to the equations of motion in terms of potentials. 



In the fourth section, we observe that the Lagrangian density can be reduced 
to some expression of the trace of the stress tensor. 

We extend the results of Seliger and Whitham by writing that / tr a dxdt 

is stationary, where tr a stands for the trace of the stress tensor. Now, we 
consider tr a as a function of the enthalpy h, the entropy s, the material 
variable X and of the tensor C where C is the right Cauchy-Green tensor. 

The variations are taken with respect to variables s, X , 0, ip, 3 introduced in 
the same way as in relation (5) and under the constraint: 

h - \~V 2 + </> + #- X T H + n = 0. 

Three appendices present auxiliary calculations and, by means of the connec- 
tive derivation associated with the velocity field, show the equivalence between 
Euler equations and our new motion equations in terms of potentials for isen- 
tropic processes. 



2 A decomposition of the stress tensor in a hyperelastic medium 



Each particle of the continuous medium is labelled by a material variable 
X, ranging over a reference configuration T>q in an Euclidian space [10]. The 
reference density p is given as a function on V [11]. 

The expression x = 0(X, t) of the spatial position describes the motion of the 
continuous medium. Generally, <f>(., t) is a twice continuously differentiable 
diffeomorphism of V onto a compact orientable manifold T> t constituting the 
image of the material at time t (see Appendix 1). As usual, we denote by F 
the deformation gradient; then, C = F T F is the right Cauchy-Green tensor. 
Recall that: 

p(detC)^=p (X). (6) 

The internal energy density is supposed to be a function of the tensor C, the 
specific entropy s and the material variable X (the hyperelastic medium is not 
necessary materially homogeneous): 

e = e(C,s,X). 

In the case of an isentropic process, this leads us back to a medium constituted 
with hyperelastic points (see Ref.10, p. 19). Now, let us put: 

1 

Q> — (J 

(detC)3 



Then, C = — C and the internal energy density can be expressed in 



P 
the form: 



e = f(p,C',s,X). (7) 



When / is independent of both C and X, we retrieve the case of an elastic 
materially homogeneous fluid. 

Let us observe that the independent variables p and C are substituted to C. 
Since essentially det C = 1, the variable p corresponds to the change of volume 
while the tensorial variable C represents the distorsion of the medium. This 
point is fundamental for the decomposition of the stress tensor and will be the 
key of the representation of the equations of motion in terms of potentials. But 
/ is defined on the manifold det C — 1 and it is more convenient to introduce 
the function g such that 

e = g(p, C, 8 , X) = f(p, f t , s, X) (8) 

(det C)3 

where g is a homogeneous function of degree zero with respect to C. We assume 
g to be a differentiable function of p, C, s, X. The stress tensor can be written 
[11,12]: 

de 

a = 2pF dc pT - (9) 

From (6) and (8), it follows: 

Using Eq. (27) proved in Appendix 2, we deduce: 



a 



2 dg n _„dg 



° = -"i I+2 " F ^ FT ' < 10) 

where / represents the identity tensor. Hence, with (10) and (29), we deduce: 

a = —pi + o\ 



with: 



2 dg _ n _„dg 



P = p 2 jf, a 1 = 2pF-^F T and tra ± = 0. (11) 

op oC 



3 A new transformation of the Hamilton principle for isentropic 
processes in elastodynamics 



The Lagrangian density in the classical principle of Hamilton has the form [2]: 

We introduce the variables C and s by means of (8) and the various quantities 
p, V, C, s and Q involved in the Lagrangian density become functions of X 
and of the motion <p of the continuous medium [2] . 

Let us denote by 5, the virtual displacement as defined by Serrin (see Ref. 2, 
p. 145). Under the Eulerian form, the variational principle of Hamilton reads: 

For all variations 5, vanishing on the boundary of W, one has 

5 L dxdt = 0. 
Jw 

Consequently, equations of the motion of isentropic processes in elastodynam- 
ics are obtained. 

In another way, it is convenient to consider the quantities p,V,s and X as 
variables. To be able to take into account the constraints imposed on p, s 
and X, it is necessary to introduce appropriate Lagrange multipliers [2,7,8] (C 
being a function of X by its derivative with respect to X). We are led to the 
variational principle: 

Theorem 1 Let <fi, if) and S be three Lagrange multipliers where <fi, ip are 
scalars and 3 is a vector, the conditions for 

I I P(-V 2 -e-Q) + <f>{-£- + div pV) - pijs - pH T X 1 dxdt = (12) 
Jw [2 at J 

to be stationary for every variation of p, *V,s,X, </>, ip, S vanishing on the bound- 
ary ofW yield the equations of isentropic motion and the constraints (2)-(4)- 

The Lagrange multiplier H is associated with Lin's constraint: 

^( P X^) + div(pVXO = 0. 

The latter arises from the fact that the Lagrangian coordinates no longer need 
to be given by explicit expressions; only the velocity field V must be such that 
it should be possible to obtain such a coordinate system by integration [5,7,8]. 



By integrating by parts the two expressions in integral (12) 

6 ( — — h divpV and — pS T X, 
\dt J 

with zero- variation for the terms on the boundary of W, we are reduced to an 
equivalent theorem: 

Theorem 2 With p, V, s, X, 0, ^, H as independent unknown functions with 
variations vanishing on the boundary ofW, and 

A = p(Iv 2 -e-fi-0-^s + X r H) (13) 

£/ie conditions for 

A da: eft (14) 



to fre stationary yield the equations of isentropic motion in elastodynamics and 
the constraints (2)-(4)- 

Proof of the theorem: 

d 
In the last part of this paper, — — will be used to denote the gradient associated 

ox 
with a tensor quantity defined on the motion space V t and ( )' will be used 

to denote the partial derivatives. Let us write 9 — g' s , the Kelvin temperature, 

and 

h = e + pg' p . 

We call h, the scalar specific enthalpy density and shortly enthalpy in this 
paper. The variations corresponding to V, p, s and X in expression (14) lead 
to the following equations: (see Appendix 3) 

«V: v = £ + «£-3r£ (15) 

OX ox ox 

Sp: A^ = 0, hence: ^V 2 - h - Q - - i/js + X r H = (15 2 ) 

5s: -p6 + pij = (15 3 ) 

oX: S = g'i + -[div(a 1 F)Y. (15 4 ) 

Then, the variations of <p, ip and H give the constraints: 

-£ + divpV = 0, s = and X = 0. 
at 

By denoting A = <fi — H r X, we deduce the equations in terms of potentials for 
the motion: 

V* = £ + <%■ + STF-' 

OX ox 



A = iv 2 -/i-fi, ^ + divpV = 0, 

2 dt 

^ = 0, s = 0, (16) 

S r = <?x + - [div(ffiF)], X = 0. 

P 

In Appendix 1, we prove directly the equivalence between system (16) and the 
Euler equations joined with constraint requirements. 



4 A simplified form of variational principle 



In the same way as Selinger and whitham, in the case of perfect compressible 
fluid [8], we shall show that the form of the variational principle given in 
Section 3 may be considerably simplified. We obtain a similar result in the 
case of isentropic motions for elastodynamics with large deformations. 

Equation (13) may be written A = ph' + p 2 g' hence A = ph.' + p, where 
p — — tr a. From (15 2 ), we obtain for the actual motions of the medium: 

h=p. 

The Legendre transformation of pg with respect to p gives — p with p being a 
function of the variables h, s, C, X, 

p = l(h,s,C,X), 

where I is a homogeneous function of degree zero with respect to C. 

Let us examine a converse. First, we notice that, for the thermo-elastic points 
of the medium, the knowledge of p as a function of the enthalpy h, the en- 
tropy s, C and X, homogeneous of degree zero with respect to C, permits the 
deduction by means of conventional thermodynamic relations, of the values of 
mass density, temperature and internal energy density. Then, 



and 



By differentiating (18), we obtain 
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= Pl 




9 = - 


Ph 




e = 


h- 


p_ 

Ph 
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;ain, 








ds 


= dh- 


dp 

P 


P 2 


dp. 



(17) 



Or 

1,1 1 1 p 

de = dh p' h dh p' ds p'-v rfX tr(p c dC) H — - dp. 

P ' P P P P 2 

From (17), we deduce: 

P 1 1 

de = 9 ds H — ^ dp p'xdX. tr(p r dC). 

P P P 

Recall that the internal energy density £ is in the form: 

£ = g(p,C,s,X), 

where g is a homogeneous function of degree zero with respect to C. If we 
choose dX.— 0, ds = and dC = Cdr, from relations tr(p' c C) = 0, tr(g' c C) = 
0, and (28) we obtain: 

P 2 #p = P 
(see Appendix 2). 

We can write the thermodynamic relations 

p 2 g' P = p, pg'c = -p'c, d = 9s, pg'-x = -p'x- 

The velocity vector may be written under a Clebsch representation: 

dx dx dx ' 

where <p, if) are scalar and H is a vector, constituting evidently a redundant set 
of unknown functions. 

dp 

As stated in the last principle, where the equation — — h div pV = specifies 

ot 

p, s = specifies s and X = specifies X, we shall require h to satisfy 

h = -V 2 -i-^s + X T S-0. (20) 

(This constraint is connected with partial result (15 2 )). We may write: 



Theorem 3 For every variation of s,X, ip, <f), 3 submitted to 
h - -V 2 + '<j> + ij) s - X r H + n = 

and vanishing on the boundary ofW, the conditions for / tr a dxdt to be 
stationary yield the equations of isentropic motions in elastodynamics and the 



relations (2)- (4) (the variations of C are deduced from the variations o/X by 
means of (30)). In the same way as in paragraph 3, we obtain System (21): 

50: ^ + divpV = 0, (21) 

Sip : s = 0, 

5s : $ = 9, 

5B: X = 0, 

SX: H r = #x + - div(o-iF). 

P 

TTie sei -Egs. (19-20) and System (21) constitutes a system of equations of the 
motion and side conditions. 

In this new form of variational principle, constraints are different from those 
of Hamilton's principle. The result is that Lagrangian and independent pa- 
rameters are completely changed. 

Let us notice that with the new parameter A = (f> — H T X and with a Clebsch 
representation for the velocity vector in the form 

ox ox 

we may write: 



Theorem 4 For every variation of s,X,ip, A, H submitted to 

h - -V 2 + A + i\) s + S r X + Q, = 

and vanishing on the boundary of W , the conditions for I tr a dxdt to be 

Jw 
stationary yield the equations of isentropic motions and the relations (2)-(4). 



5 Conclusion 



The variational statement in Theorem 4 is far from Hamilton's principle . The 
resulting motion equations involve thermodynamic variables like temperature 
and entropy. 

The constraint condition h — — V + A + ip s + S T X + Q = involves the 

specific enthalpy. Although this constraint explicitly uses entropy and the non 
uniquely defined parameters A and if), the variational principle in this form 



10 



is interesting in view of the fact that the stress tensor a has an experimental 
meaning. 
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6 Appendix 1 

Equations of isentropic motions in elastodynamics 

Preliminaries [13,14] 
The motion of the medium consists in the ^-dependent C 2 -diffeomorphism 



X 



X 1 
X 2 
X 3 



eVn 



X 



X 



X 



XX 



eV t , [11,12]. 



We immediately deduce F = -7—F and F~ l = —F^ 1 ——. Let us write T*(V t ) 

ox ox 

for the cotangent fiber bundle of V t and T*(V t ) the cotangent linear space to 

T> t at x; then, 

xev t — ► L(t,x) eT*(v t ) 

represents a differential form field on T> t . 

Let us write T*(D ) the cotangent fiber bundle of V and T^_(V ) the cotangent 
linear space to T>q at X. The mapping <pl is induced by <j) t for the form fields. 
The convective derivation d c of a form field L is deduced from the diagram: 



L e T*(V t 

14 



LF e T*(V ) 

I A 

* dt 



(22) 



L + L— G T*(V t ) £- LF + L— F e T*(V ) 
ox ox 



dV 
where L + L— — is the Lie derivative of L with respect to the velocity field V 

ox 
that is the infinitesimal transformation of the one-parameter group of trans- 
formations (f) t [15]. 

Consequences 

Let b be a scalar field on T> t assumed to be an Euclidian space; grad represents 
the gradient operator on T> t . We define two form fields by their values V r and 
(grad6) r . From (22) we deduce: 



4(V^ + A ( V) 



(23) 



12 



(where T is the acceleration), and 

4(grad6) T = (gradfe) T (24) 

d^LF- 1 ) = LF-\ (25) 

Potential equations 
With the notations of Section 2, the motion equation is: 



^ + g-div„ + ^0. (26) 



Relation (8) gives: 



de = 6 ds + -^ dp + g' x dX + tr(g' c dC) . 

By symmetry property of the tensor g' c , we obtain: 

tr(g' c dC) = 2tr(g' c F T dF), 
and Eq. (11) implies 

tr{g' c dC) = -trtaxdFF- 1 ). 
p 

From the definition of the specific enthalpy, we deduce: 

1 p 

grad h = grad e H — gradp grad p, 

p p 2 

f = e* + I* + A *-> + Idiv^FjF- - Idiv,, 

ox ox pox p p 

Equation (26) becomes: 

pT^ + p-^-(h + n)-div(a 1 F)F- 1 - p6^- g' x F- 1 = 0. 

Equation (23) implies: 

c) ~\ (9s 1 

4(Vr) = ^ ( 2 ^ ~ k ~ n) + ^ + " ^ div ( aiF ) + P&} p ~ 1 - 

Let H T be a form field such that H T = — div(<7iF) + g'-^. With substituting 
from Eq. (25), 
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Let A and ip be two scalar fields verifying A = -V 2 — h — Vt and ip = 6. By 
adding to S T , if necessary, an appropriate one- form with a zero convective 
derivation, Eq. (24) and s = (isentropic process), give: 

OX ox 
where 

A = -V 2 -h-n, i) = Q and S = - div(aiF) + g'^. 
1 p 

These equations along with the relations (2), (3), (4) are potential equations 
of the motion such as have been obtained in paragraph 3. 



7 Appendix 2 

Some calculus 

Mass density 
By differentiating Eq. (6) and using Jacobi's identity, we obtain: 



Hence, 



dp = - p tr(C- L dC) + (det C)~ t^- dX. 
2 oJ\- 



P 'c = -\pC- 1 . (27) 



If we choose dC = Cdr and dX.— 0, this reduces to 

dp = ~pdr. (28) 

With respect with C, g is a homogeneous function of degree zero. We deduce 
immediately (Euler identity): 

tr(g' c C) = 0. (29) 



Cauchy- Green tensor 
Let us calculate the Cauchy-Green tensor variations: 

dC = dF T F + F T dF, 5F = -F^—F. (30) 

ox 



14 



Hence, 



sc = - 



c^ F+( A r 



dx 



dx 



(31) 



8 Appendix 3 

Proof of Eqs. (15) 

Let us write (14) under the form: 



A<fa tt=/Jiv'- e -n-g-£v-4 

w Jw [2 at ox at 

3s <93 <9S ) 

-^ v + xr ar +x '* v \ A * dt 



For every £V-variation, we deduce: 

ds ,„<9E 



Jw { 



ijj-+X T —\ 6Vdxdt = 0, 
ox ox ox 



and we obtain immediately Eq. (15). An integration by parts, on the boundary 
of W, leads with the same calculation to Eqs. (15 2 ),(15 3 ) and the constraint 
relations. 

For every (DC-variation, vanishing on the boundary of W, we obtain: 

-pg' x 5X + pS 5X- ptr(g' c 5C)\ dxdt = 0, 

where SC is a function of 5X (see Appendix 2, Eq.(30)). 
By using Eqs. (11), (30), we obtain: 

,<9<5X\ 



w 



-ptr{g' c 5C) = tr[a 1 F- 



dx 



Moreover, we can write: 



Hence 



w 



w 



p(H -g^SX + tr [a x F 



d5X\ 
dx J 



dxdt = 0. 



S - t7x div (<TiF) 



5X dxdt = 0. 
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We thus obtain 

■ T 1 

s =#x + - div(o-iF). 

Let us notice that the divergence operator is applied to a tensor defined on 
£>o with values in T> t ; the partial derivatives are calculated on T> t . 

D 
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